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Abstract 

A powerful method for the description of open quantum systems is the Feshbach projection 
operator (FPO) technique. In this formalism, the whole function space is divided into two subspaces 
that are coupled with one another. One of the subspaces contains the wave functions localized in 
a certain finite region while the continuum of extended scattering wave functions is involved in the 
other subspace. The Hamilton operator of the whole system is Hermitian, that of the localized part 
is, however, non-Hermitian. This non-Hermitian Hamilton operator H e g represents the core of the 
FPO method in present-day studies. It gives a unified description of discrete and resonance states. 
Furthermore, it contains the time operator. The eigenvalues z\ and eigenfunctions (f>\ of i? e g are 
an important ingredient of the S matrix. They are energy dependent. The phases of the <f>\ are, 
generally, nonrigid. Most interesting physical effects are caused by the branch points in the complex 
plane. On the one hand, they cause the avoided level crossings that appear as level repulsion or 
widths bifurcation in approaching the branch points under different conditions. On the other hand, 
observable values are usually enhanced and accelerated in the vicinity of the branch points. In most 
cases, the theory is time asymmetric. An exception are the VT symmetric bound states in the 
continuum appearing in space symmetric systems due to the avoided level crossing phenomenon 
in the complex plane. In the paper, the peculiarities of the FPO method are considered and three 
typical phenomena are sketched: (i) the unified description of decay and scattering processes, (ii) 
the appearance of bound states in the continuum and (iii) the spectroscopic reordering processes 
characteristic of the regime with overlapping resonances. 
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I. INTRODUCTION 



An exact description of open quantum systems meets the problem to consider simulta- 
neously the wave functions of discrete and scattering states. Both types of wave functions 
are completely different from one another. The discrete states A characterize the spectrum 
of the system and are normalized according to the Kronecker delta 5\\> while the scatter- 
ing states are continuous in energy E and can be normalized according to the Dirac delta 
function 5(E — E'). The wave functions of discrete and scattering states appear in a com- 
bined manner in most physical expressions characteristic of open quantum systems. Special 
mathematical considerations are necessary therefore in order to receive physical values. 

In the iV-level Friedrichs model , the total Hamiltonian H is defined by 

H = H + pV (1) 

where p is a real number and Hq is the so-called free Hamiltonian 

H = 'S2u n \n){n\ + / uo\uo) {uo\ p{uj)duo . (2) 

Here, \n) and \uo) satisfy the orthonormality condition (n\n f ) = 5 nn > and (uj\uj') = 5{uo — 
u')/p(uj), and (n\ou) = 0. The sum runs over the (finite) number of discrete basic states 
\n) and the integral is over the considered energy region with K u = {uj\p{oj) ^ 0}. The 
interaction Hamilton operator V describes the coupling between \n) and \u), 

V = J2 {y n {uj) \u)(n\ + v* n (u) \n)(uyj p(uj)du (3) 

where v n (uj) is the interaction matrix element between \n) and \uo). In the Friedrichs model, 
the Schrodinger equation with the Hamiltonian H is directly solved. It is however not easy 
to receive results that are of physical interest in a broad range of parameters. For an example 
of the troubles see the study on bound states in the continuum Q]. 

Another method to solve the Schrodinger equation with the Hamilton operator ([1]) is 
the use of the Feshbach projection operator (FPO) technique. In this method, the basic 
equations for the wave functions of the states A and to are solved separately such that the 
main problem of the Friedrichs model is avoided. In the FPO formalism Jjjj], the full function 
space is divided into two subspaces: the Q subspace contains all wave functions that are 
localized inside the system and vanish outside of it while the wave functions of the P subspace 



are extended up to infinity and vanish inside the system, see [l]. The wave functions of the 
two subspaces can be obtained by standard methods: the Q subspace is described by the 
Hermitian Hamilton operator Hb that characterizes the closed system with discrete states, 
while the P subspace is described by the Hermitian Hamilton operator He that contains the 
continuum of scattering wave functions. Thus, H = Hb + He in (121). The coupling matrix 
elements are calculated according to (j3J) by using the eigenf unctions A of Hb instead of the 
basic wave functions n that appear in ([3]). An example for the difference between the A and 
the n is the following: in nuclear structure calculations, the n are determined by the Slater 
determinants while the A are the shell-model wave functions. Furthermore, // = 1 in (pQ), 
i.e. there is no free parameter in H in the framework of the FPO formalism. All values are 
determined by fixing the potential and the coupling matrix elements between the different 
discrete and scattering states and by defining the two subspaces in a meaningful manner. 

In the FPO formalism, the closed system (defined by the Hamilton operator Hb) will be 
opened by coupling the wave functions of the Q subspace to those of the P subspace under 
the assumption P + Q = 1. Due to this coupling, the discrete states of the closed system 
that lie above particle decay thresholds, pass into resonance states of the open system. The 
states below decay thresholds receive, as a rule, some energy shift but remain discrete. The 
resonance states have, in general, a finite life time. 

The FPO method is introduced by Feshbach [3] forty years ago in order to describe nuclear 
reactions. At that time it was impossible to perform all the calculations in the two subspaces 
as well as those for the coupling matrix V. Instead, Feshbach used statistical assumptions 
for the narrow states of the Q subspace (compound nucleus states) and treated exactly only 
the so-called direct (fast) reaction part. In this manner it was possible to formulate a unified 
description of nuclear reactions, i.e. of the fast direct nuclear reaction part and the much 
slower compound nuclear reaction part. 

In the present-day calculations on the basis of the FPO method, all the calculations in the 
Q subspace are performed with the same accuracy as the calculations for the corresponding 
closed system (Q = 1). Also the coupling matrix elements v\(w) are calculated. These 
calculations represent therefore a unified description of structure and reaction phenomena 
{4]. They allow to draw general conclusions on the behavior of open quantum systems under 
different conditions, i.e. by controlling them in a broad parameter range. 

In the present paper, the FPO method is sketched (Sect. 2) and some basic peculiarities 
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are discussed (Sect. 3). An important feature of the method is the appearance of the non- 
Hermitian Hamilton operator H c g in an intermediate stage of the model. This Hamilton 
operator allows a unified description of discrete and resonance states. The complex eigen- 
values z\ and eigenfunctions (f)\ of H c q determine decisively the S matrix (being unitary for 
all parameter values) and observables as will be shown by means of concrete examples in 
the following sections. Typical features of the formalism are time asymmetry (Sect. 4), the 
existence of bound (T symmetric) states in the continuum when the system is V symmet- 
ric (Sect. 5) and spectroscopic reordering processes that take place under the influence of 
branch points in the complex energy plane in the cross over from the weak-coupling regime 
to the strong-coupling one (Sect. 6). 

II. BASIC RELATIONS OF THE FESHBACH PROJECTION OPERATOR (FPO) 
FORMALISM 

The basic equation of the FPO formalism 

(H-E)^ = (4) 

is solved in the whole function space by dividing it into the two subspaces P and Q with 
P + Q = 1. The Hamilton operator H is Hermitian. It contains the decay of the subsystem 
localized in the Q subspace, into the surrounding P subspace where the decay products can 
be detected. The excitation of the states localized in the Q subspace may take place via one 
of the channels C included in (T4j) or by another process that can be described by a source 
term F appearing on the right-hand side of (jl]), for details see Q]. 

In solving (j4]) in the whole function space P + Q = 1 by using the FPO technique, the 
non-hermitian Hamilton operator 

H eS = H B + J2V BC „ + 1 „ V CB (5) 

c h,^ — ti c 

appears which contains Hb as well as an additional non-hermitian term that describes the 
coupling of the resonance states via the common environment. Here Vbc, Vcb stand for 
the coupling matrix elements between the eigenstates of Hb and the environment [4| that 
may consist of different continua C. The operator H e g characterizes the part of the problem 
that is localized in the Q subspace while the operator H describes the problem in the whole 
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function space P + Q. Therefore, H e g is non-Hermit ian and H is Hermitian. The operator 
H e Q is explicitly energy dependent and symmetric, 



{H eS -z x )(P x = 0, (6) 

its eigenvalues z\ and eigenf unctions <ft\ are complex. The eigenvalues provide not only the 
energies of the resonance states but also their widths. The eigenfunctions are biorthogonal. 
For details see [4j. 

The eigenvalues and eigenfunctions of Hb contain the interaction u of the discrete states 
which is given by the nondiagonal matrix elements of Hb- This interaction is of standard 
type in closed systems and may be called therefore internal interaction. The eigenvalues and 
eigenfunctions of H c g contain additionally the interaction v of the resonance states via the 
common continuum (v is used here instead of the concrete matrix elements of the second 
term of H c g). This part of interaction is, formally, of second order and may be called external 
interaction. While u and Re(w) cause level repulsion in energy, Im(t> ) is responsible for the 
bifurcation of the widths of the resonance states (resonance trapping). The phenomenon 
of widths bifurcation (resonance trapping) has been proven experimentally in microwave 
cavities p. 

Since the effective Hamilton operator §5$) depends explicitly on energy E, so do its eigen- 
values z\ and eigenfunctions 0a- Far from thresholds, the energy dependence is weak, as a 
rule, in an energy interval of the order of magnitude of the width of the resonance state. 
The solutions of the fixed-point equations E\ = Re(z\)\E=E x an d of T A = —2Im(z\)\ E=Ex 
are numbers that coincide (approximately) with the poles of the S matrix. In the FPO 
formalism, however, it is not necessary to consider the poles of the S matrix since the spec- 
troscopic information on the system follows directly from the complex eigenvalues z\ and 
eigenfunctions <p\ of H e g. Moreover, in the physical observables related to the S matrix the 
eigenvalues z\ with their full energy dependence are involved, see (|T5|) . Due to this fact, in- 
formation on the vicinity (in energy) of the considered resonance states such as the position 
of decay thresholds and of neighboring resonance states is involved in the S matrix and can 
be received. Such an information can not be obtained from the poles of the S matrix being 
(energy-independent) numbers. 

In contrast to the trajectories Z\(X) of the eigenvalues of a Hermitian Hamilton operator 
(where X is a certain parameter), those of a non-Hermitian one may cross in the complex 



plane. The crossing points are branch points (called exceptional points in the mathematical 
literature). Physically, they are responsible for the avoided level crossing phenomenon ap- 
pearing in their vicinity. More precisely: in approaching the branch points under different 
conditions, we have level repulsion (together with widths equilibration) or widths bifurcation 
(together with level attraction), see |4(. 

The eigenfunctions <p\ of H e R are complex and biorthogonal. The normalization condition 
(0 A ft |0 A ght ) = (4>*\\4>\) fixes only two of the four free parameters 7]. This freedom can 
be used in order to provide a smooth transition from an open quantum system (with, in 
general, nonvanishing decay widths T\ of its states and biorthogonal wave functions 0a) to 
the corresponding closed one (with r A — > and real wave functions that are normalized in 
the standard manner): (4>\\4>\) — > (4>\ \ <px) = 1 if the coupling vectors in the non-Hermitian 
part of (jSJ) vanish. That means, the orthonormality conditions can be chosen as 

(0aI<M=V' (7) 

with the consequence that [4j 

(0 A |0 A > = A\ > 1 (8) 
B* = (0 A |0 AYA ) = -B w ee - (0 AYA |0 A ) ; \B%\ > . (9) 

Approaching the branch point where the two eigenvalues z\ and zy coalesce, A\ — > oo and 
|.B A '| — > oo. The normalization condition ([7]) entails that the phases of the eigenfunctions 
in the overlapping regime are not rigid: the normalization condition (0 A |0 A ) = 1 is fulfilled, 
in this regime, only when Im(0 A |0 A ) oc Re <p\- Im A = 0, i.e. by rotating the wave function 
at a certain angle (3\. For details see js}. The phase rigidity defined by 
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(0a|0a) (Re0 A ) 2 + (Im0 A ) 2 A\ 1 ] 



is a useful measure [9j for the rotation angle (3\. When the resonance states are distant from 
one another, it is r A w 1 due to (0 A |0 A ) ~ (0 A |0 A ). In approaching a branch point in the 
complex energy plane 0, [h]], we have r A — > 0. Therefore 1 > r A > 0. 

The phase rigidity r A is a measure for the degree of alignment of one of the overlapping 
resonance states with one of the scattering states £q of the environment. This alignment 
takes place at the cost of the other states that decouple, to a certain extent, from the 
environment (widths bifurcation or resonance trapping [4]). It agrees with experimental data 



[ill ], according to which the phase rigidity drops smoothly from its maximum value r\ — 1 
far from the branch point to its minimum value r\ = at the branch point, see 



The solution of (HI) reads 



4] 



l^) = l^) + EI^) i f^) (ii) 

where 

\£l c x ) = (l + E+ l _ H V CB )\<Px) (12) 

is the wave function of the resonance state A and the £j£ are the (coupled) scattering wave 
functions of the continuum into which the system is embedded. The expression ( TTT1) follows 
by applying P + Q = 1 to (j4j) without any approximations. It is therefore an exact solution 
of (jl]). The representation of (being solution of (jlj) with the Hermitian operator H) 
in the set of wave functions (f)\ (being solutions of ([6]) with the non-Hermitian operator 
H c h) characterizes the consideration of localized states in the FPO formalism. According 
to ffTTj) . the eigenfunctions <p\ of the non-Hermitian Hamilton operator H e g give the main 
contribution to the wave function in the interior of the system, 

l*g>-l*S> = £<&lfc>; c£ A = if^ (13) 

A ^ 

and 

- = £<& (0f I = £<& (011 • (14) 

A A 

The weight factors c^ A and d^ x (= in the scattering process) contain decay and excita- 
tion, respectively, of the states A at the energy E. 



The S matrix is unitary. It can be obtained from (^\V\^/q), see [4|. The amplitude of 
the resonance part we are interested in, is given by 

^ = iZ<£\VM& = <Z ® lV M™ Vl ® . (15) 

A A Z\ 

The resonance structure of S ves is determined by the eigenvalues z\ of H e g as long as the 
resonances do not overlap (see Sect. 6). 
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III. PECULIARITIES OF THE FESHBACH PROJECTION OPERATOR (FPO) 
FORMALISM. UNIFIED DESCRIPTION OF DISCRETE AND RESONANCE 
STATES 



The characteristic features of the FPO formalism consist, above all, in the fact that the 
solution in the whole function space can be represented in the set of wave functions {0a} 
that describe the localized part of the problem. The localized wave functions represent a 
subspace of the whole function space with the consequence that the corresponding Hamilton 
operator H e g is non-Hermitian. 

The main advantages of the FPO formalism consist in the following. 

(i) The spectroscopic information on the resonance states is obtained directly from the 
complex eigenvalues z\ and eigenfunctions <p\ of the non-Hermitian Hamilton operator H c r. 
The z\ and <p\ are energy dependent functions, generally, and contain the influence of neigh- 
boring resonance states as well as of decay thresholds onto the considered state A. This 
energy dependence allows to describe decay and resonance phenomena also in the very 
neighborhood of decay thresholds and in the regime of overlapping resonances. Since also 
the coupling coefficients between system and continuum depend on energy, the unitarity of 
the S matrix is guaranteed for all parameter values, see e.g. 



(ii) The resonance states are directly related to the discrete states of a closed system 
described by standard quantum mechanics (with the Hermitian Hamilton operator Hb)- 
They are generated by opening the system, i.e. by coupling the discrete states to the 
environment of scattering states by means of the second term of the Hamilton operator H e ^. 
Therefore, they are realistic localized (long-lived many-particle) states of an open quantum 
system. 

(iii) The properties of branch points and their vicinity can be studied relatively easy. 
At these points, two (or more) eigenvalues z\ of H c g coalesce. Since it is not necessary to 
consider the poles of the S matrix in the FPO formalism, additional mathematical problems 
at and in the vicinity of branch points (exceptional points) in the complex plane are avoided. 

(iv) The phases of the eigenfunctions <fi\ of H c r are not rigid in the vicinity of a branch 
point. This fact allows to describe the spectroscopic reordering processes in the system that 
take place under the influence of the scattering wave functions of the environment into 
which the system is embedded. 




S 



These features are involved in all present-day [13( calculations performed on the basis of 
the FPO formalism. In numerical studies, the main problem arises from the definition of 
the two subspaces Q and P such that it is meaningful for spectroscopic studies (the cross 
section is, of course, independent of the manner the two subspaces are defined, for more 
details see ^]). The basic idea is the following: H B describes the closed system (localized 
in the interior of the system) which becomes open when embedded in the environment of 
scattering wave functions described by Ho- Therefore, all values characteristic of resonance 
states can be traced back to the corresponding values of discrete states by controlling the 
coupling to the continuum. That means with v — > 0, the transition from resonance states 
(described by the non-Hermitian H ef[ ) to discrete states (described by the Hermitian H B ) 
can be controlled. 

As to the mathematical properties of branch points in the complex energy plane, it will 
be mentioned here only that the phase jump behavior of the wave functions at and in the 
vicinity of the branch points and their topological structure are investigated in the framework 
of the FPO formalism by means of numerical calculations for some special cases as well as 

□ n n n 

in analytical studies [10, ll5|, ll6j • The topological structure of the branch (exceptional) 



points in the continuum differs from that of diabolic points. The results obtained in the FPO 
formalism are in full agreement with those of rigorous mathematical studies for a symmetric 
non-Hermitian 2x2 matrix Hamiltonian At the exceptional (branch) point, the two 
different right and left eigenvectors of the non-Hermitian Hamilton operator are linearly 
dependent, \(f>\) <-> ± i\(f>\>); (</>x\ <-» T ^(0a'I> an d are supplemented by the corresponding 
associated vectors defined by Jordan chain relations. Moreover, the results are in agreement 
with those of experimental studies on microwave cavities 1171, see Ref. ffl. The cross section 

n n 

(S matrix) behaves smoothly at the branch point [4|, [18|. 

The physical meaning of the branch points in the complex energy plane is based upon 
their topological structure and their relation to the phenomenon of avoided level crossing, 
i.e. to, respectively, level repulsion and widths bifurcation occurring in approaching them 
under different conditions. Level repulsion is accompanied by widths equilibration, while 
widths bifurcation occurs together with level attraction (formation of clusters). For further 
details see e.g. j^]. 

In the FPO formalism, a unified description of discrete and resonance states is involved. 
First, the resonance states pass smoothly into discrete states by reducing the coupling 
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strength v between system and environment. The smooth transition is guaranteed by choos- 
ing the orthonormality conditions for the eigenfunctions <p\ of H e s according to (j7]). Sec- 
ondly, the localized states in the Q subspace are eigenstates of H eS . Generally, they may 
lie above as well as below the particle decay thresholds. In the first case, their widths T\ 
are different from zero (with the exception of bound states in the continuum, see Sect. 5) 
while they are zero in the second case. In the last case, H e g is real: the residuum of the 
second term of (jSJ) vanishes and only the principal value integral is different from zero. It 
causes, generally, some shift of the position of the state A, A(E^ — E\) ^ 0, see {4]. Also 
the wave functions of these states differ, generally, from those of the eigenfunctions of Hb, 
<fi\ 7^ That means, discrete and resonance states are described in a unified manner 
by the non-Hermitian Hamilton operator H e ^\ not only the resonance states lying above 
particle decay thresholds, but also the discrete states below the thresholds are influenced by 
the continuum. As a consequence, the branch points in the complex plane determine also 
the properties of bound states. The most impressive example is the phenomenon of avoided 
crossing of discrete states known since the very beginning of quantum physics [19[. It can 
be traced from discrete states up to the branch points in the complex plane by parameter 



variation 



15] . At the branch points, nonlinearities play a role. The branch points intro- 



duce therefore nonlinear effects into quantum mechanics 4J. This fact may lead to a deeper 
understanding of the relation between avoided level crossings and quantum chaos. 

These results show that the non-Hermitian Hamilton operator H e s is basic for the de- 
scription of the localized states in realistic quantum systems. It describes the spectroscopic 
properties of discrete and resonance states in a unified manner. The restriction to the 
Hermitian Hamilton operator Hb [see in the standard quantum mechanics is an ap- 
proximation: the principal value integral of the second term of (j3J) is effectively taken into 
account, however the residuum is neglected. Although this approximation works well in 
very many applications, it leaves open some fundamental questions. In these cases, the 
non-Hermiticity of the Hamilton operator H e g (and the nonlinear effects related to it) can 
not be neglected. 

Another peculiarity of the FPO formalism is the existence of a time operator which is the 
residuum of the non-Hermitian Hamilton operator H e ^. The life time T\ of a resonance state 
follows from the eigenvalue z\ of H c q in the same manner as the energy E\ of this state. 
Both values are fundamentally different from the time t and the energy E. They characterize 
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the state A while t and E appear as general parameters. In the closed system with the 
Hermitian Hamilton operator Hb, only the energies Eb of the states can be determined. 
The eigenvalues are real and the widths are zero, Tb = 0. Due to the coupling to the 
continuum, energy shifts E x — Eb of the states appear as well as the finite life times T\ oc 
(T\ — r^) _1 = T^ 1 of the resonance states. Both, the energy shifts and the finite life times, 
follow from the non-Hermitian coupling term of if e ff [see Eq. Usually, the numbers Ex 
and Tx can be obtained directly from the z\. Only in the case the zx are strongly dependent 
on energy, the corresponding fixed-point equations have to be solved. The energies Ex 
and life times Tx of the resonance states A of an open quantum system are bounded from 
below (see jg] for the discussion of the brachistochrone problem in open quantum systems). 
Mathematically, the existence of the time operator entails the time asymmetry involved in 
the FPO formalism. 



IV. UNIFIED DESCRIPTION OF RESONANCE AND DECAY PHENOMENA 

The time dependent Schrodinger equation reads 

H efi if E (t) =ih^ * E (t) . (16) 

The right solutions may be represented, according to ffl3|) . by an ensemble of resonance 
states A that describes the decay of the localized part of the system at the energy E, 

(right) ^ __ e -iH cB t/h (right) ff^^ 

= E e-^ h cxo \4 igkt) ) (17) 

A 

with |0^ lght ' ) ) = \4>x) an d cxo = (<t>Vy\£>c) I — -2a)- The zx and 4>x are the (energy depen- 
dent) eigenvalues and eigenfunctions of the time- independent Hamilton operator if eff , Eq. 
([5]), while the are the scattering wave functions of the environment. The left solution of 
(flEj) reads 

= E i^ldxte^ (18) 

A 

with (0 A left) | = (<f)*x\ and d X t = c* xo = (£<f| V\<f>x)/(E - z* x ) for the scattering process. It 
describes the excitation of the system at the energy E. For other excitation processes, 
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e.g. via a source term F on the right-hand side of PJ, see j^]. Here d\ t is, generally, time 
dependent. 

By means of (TTTf) and ffTgj) the population probability 

(*g (left) (*)|^ (risht) (^)) = E c ao e" r ^ (19) 

A 

at the energy E can be defined. The decay rate reads 

U*) = -|ln(^ (left) (*)|^ (risht) (t)> 
1 EaTa c A0 d Af e~ r ^ h 

For an isolated resonance state A, (1201) passes into the standard expression 

fcgr(t) -> k x = T x /h. (21) 

In this case, the value k\ is constant in time and corresponds to the standard relation 
t\ = h/Y\ with r\ = l/k\. It describes the idealized case of an exponential decay law and, 
according to (1151) . a Breit-Wigner resonance in the cross section. Generally, deviations from 
the exponential decay law and from the Breit-Wigner line shape appear under the influence 
of neighboring resonance states and (or) of decay thresholds. Also the background term 
appearing in most reactions may cause deviations from the ideal exponential decay law. For 
details see ||| 

The expressions (I17j) and ( 1T81) are valid only when (fl3l) holds, i.e. at times t at which 
the wave functions have a localized part in the interior of the system at the energy E so 
that the representation (1131) is meaningful at this energy. According to (1171) and (I18p . this 
is the case for times t > to where to is a finite value. Without loss of generality, it can be 
chosen t = 0. The quantum system described in the framework of the FPO formalism is 
therefore time asymmetric. The time asymmetry is involved in the non-Hermitian part of 
the Hamilton operator H c r (which contains the time operator), as can be seen immediately 
from the expression (Tl9l) for the population probability. 

The consideration of only the time interval < t < oo in ffl6l) is related to the fact that 
the decay of a resonance state (at the energy E of the system) starts at a finite time (say 
t = 0) at which the system can be considered to be excited, i.e. (|T3|) is meaningful at this 



energy. This fact agrees with the concept of a semigroup description introduced in 20], 
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which distinguishes between prepared and measured states. In our formalism, the decaying 
(measured) states are described by the eigenvalues and eigenfunctions of the effective non- 
Hermitian Hamilton operator H e R involved in the \^ E ), Eq. (I13p . The preparation of the 
resonance states is described by the energy-dependent (^ E \. It may be very different for 
different reactions. 

The decay properties of the resonance states can be studied best when their excitation 
takes place in a time interval that is very short as compared to the life time t\ of the resonance 
states, i.e. d\ t is a function being strongly dependent on time. In such a case, the time t = 
is well defined and no perturbation of the decay process by the stfil continuing excitation 
process will take place. In [21| . such a situation is studied in single ion experiments. The 
results demonstrate the beginning of time for a decaying state. That means, they prove the 
time asymmetry in quantum physics. 

Eq. (j20|) describes the decay rate also in the regime of overlapping resonances, see 0, [22]] . 
The overlapping and mutual influence of resonance states is maximal at the branch points 
in the complex plane where two eigenvalues z\ and z\> of the effective Hamilton operator 
ifeff coalesce. Nevertheless, the decay rate is everywhere smooth as can be seen also directly 
from (120]) . This result coincides with the general statement according to which all observable 
quantities behave smoothly at singular points. 

Another interesting problem is the saturation of the average decay rate k av in the regime 
of strongly overlapping resonances. According to the bottle-neck picture of the transition 



state theory, it starts at a certain critical value of bound-continuum coupling 23J]. This 



saturation is caused by widths bifurcation (formation of long-lived resonance states by reso- 
nance trapping by a few short-lived states [4]) occurring in the neighborhood of the branch 
points in the complex plane s|, [24J] . Widths bifurcation creates long-lived resonance states 
together with a few short-lived resonance states. The definition of an average life time 
of the resonance states is meaningful therefore only for either the long-lived states or the 
short-lived ones. The long-lived (trapped) resonance states are almost decoupled from the 
continuum of decay channels. Their widths T\ saturate with increasing bound-continuum 
coupling. The T x are almost the same for all the different states A, see [4|, i.e. T av « I\ for 
all long-lived trapped resonance states. It follows therefore 

k av « r av /h (22) 
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from (1201) . According to the average width r av , the average life time of the long-lived states 
can be defined by r av = 1/A; av . Then f )22|) is equivalent to r av = /i/T av . That means, the 
basic relation between life times and decay widths of resonance states holds not only for 
isolated resonance states [see Eq. (]2"T]) ]. but also for narrow resonance states superposed by 



a smooth background (that may originate from a few short-lived resonance states 4j, |25| ) . 
In the last case, the relation holds for the average values T av and r av . 

The expressions ([15]) for the S matrix and (1201) for the decay rate show immediately that 
the resonance phenomena (described by the S matrix) are determined by the decay properties 
of the resonance states (described by the complex eigenvalues z\ and eigenfunctions <fi\ of 
the non-Hermitian Hamilton operator H e s). Thus, the FPO formalism provides a unified 
description of resonance and decay phenomena. The expression f[T5"]) shows however also 
that, generally, the energy dependence of the eigenvalues z\ and eigenfunctions (p\ of H e g 
causes deviations from the Breit-Wigner resonance shape and the exponential decay law. 
The deviations become important for isolated resonance states due to the fact that the 



decay thresholds lie at a finite energy 26|] . They appear mainly in the long-time scale. This 



result agrees qualitatively with experimental data 271 ] . At high level density, deviations 



appear even in the short-time scale due to the mutual influence of neighbored resonance 



states, see Sect. 5 and jg]. 



V. BOUND STATES IN THE CONTINUUM (BICS) 

The question whether or not bound states in the continuum (BICs) exist in realistic 
quantum systems is of principal interest and might be as well of interest for applications. 
The reason for this interest arises from the fact that the system is stabilized at a BIC as 
well as in its vicinity, and that the wave function is localized at all times inside the system 
in spite of embedding it into the continuum of extended wave functions. 

Mathematically, the existence of bound states in the continuum is shown already in 1929 



by von Neumann and Wigner [28j. In 1985 Friedrich and Wintgen 29] considered the 
problem by using the FPO technique. They related the existence of BICs to avoided level 
crossings being another quantum mechanical phenomenon discussed by von Neumann and 
Wigner [19j in 1929. As discussed in Sects. 2 and 3, avoided level crossings are caused by 
branch points and appear in their vicinity. 
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Since BICs are states that do not decay, the population probability of these states is con- 
stant in time. This fact is called population trapping in studies on laser induced continuum 

n n 

structures in atoms [30( • Similar results are obtained [14j in the time independent approach 
by using the FPO technique and demanding a vanishing decay width for the BIC. In these 
papers, the relation between BICs and the avoided level crossing phenomenon as well as 
the stabilization of the system in a broad range of the parameter values (characteristic of 
the laser) is shown explicitly. A similar study is performed for the transmission through a 
quantum billiard where BICs appear at those energies at which the resonant transmission 



crosses a transmission zero 



31 



321 ] . Common to all these studies is the definition of a BIC 



as a resonance state with vanishing width, 

T Xo \(e=e Xo ) = . (23) 



In the FPO formalism, its energy is obtained from the solution of the fixed-point equation 
E\ = E x \ {E=Exo) . 

Generally, the relation between T\ and the coupling matrix elements is |4| 



T x = -2lm(z x ) <2rrJ2\(<Pl\V\^)\ 2 ■ (24) 

c 

This expression holds true at all energies. The inequality in ( 12^1) is caused by the biorthog- 
onality of the functions <p\ being eigenf unctions of the non-Hermitian operator §5§. As a 
consequence, a state being decoupled from all channels C of the continuum according to 

(HE\v\K) ^ Q (25) 

is a BIC with T\ = —2 Im(zA ) ~~ y [the condition ( 1251) is equivalent to (0a o I^I£c) — * due 
to the symmetry of H c g and the biorthogonality of the 0a]- The opposite case follows by 
considering the S matrix, see ffl5l) for the amplitude of its resonance part. At the position 
of a BIC, we have E — z\ — > and, due to the unitarity of the S matrix, it follows (|25|) 
for all C . That means: the decoupling from all channels of the continuum described by 
(I25p is a necessary and sufficient condition for a resonance state to be a BIC, i.e. a state 
with vanishing decay width T\ = 0. The wave function of such a BIC is, according to ( fT2|) . 
eigenfunction of H e s and, consequently, localized. 

In {33I, the advantage of the FPO method as compared to the iV-level Friedrichs model 
[2I in studying BICs for unstable multilevel systems is discussed. In contrast to the coupling 
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matrix elements (|25|) . the form factors (^cl^l^no) an d (0n o l^l£<?) considered in Q contain 
the basic wave functions 0^ of the Hamiltonian ifg of the closed system. Since the eigen- 
functions 4>\ of H c s can be represented as = J2 a \,x4>x with complex coefficients cl\ \> 
[and = J2 b\n<Pn with real b\ n and the basic wave functions 0^ of discrete states defining 
Ho according to (j2J)], a sum of individual form factors vanishes at the position of a BIC 
according to (1251) . 

£GAo,A'<£fW*>^0. (26) 
A' 

This equation is nonlinear since the coefficients cl\ \i depend on the coupling strength. It 
is therefore difficult to obtain a general solution for the multilevel case in the framework of 
the Friedrichs model. 

By using the FPO technique and demanding a vanishing decay width r^ , the multilevel 



problem can be solved. Examples are given in [14| for atoms and in [31( for quantum dots. 
At sufficiently small coupling strength v, the condition I\ = can be fulfilled only when 
the spectrum of the closed system (described by Hb) is degenerated. Concrete examples of 
BICs investigated recently analytically with the postulation T\ = for a BIC are studies 



on open quantum billiards with variable shape [32|. The relation of the BICs to the avoided 
level crossing phenomenon can be seen in all examples. The scattering phase jumps by tc in 
approaching the BIC in spite of the fact that no resonance can be seen in the cross section 

By means of the complex eigenvalues z\ of H e g, the appearance of a BIC can be traced 
as a function of a certain control parameter X, i.e. by controlling the trajectories E\(X) 
and T\(X). The BIC appears at the point X = X where T\(X ) = 0. It is even possible 
to consider the vicinity of the BIC including the cases when T\(X') is always different from 
zero and T\(Xq) corresponds to the minimum of T\(X') with a small but nonvanishing 
value T\(Xq) rs 0. This feature of the FPO technique is invaluable for applications since the 
stabilization of the system (caused by the vanishing width T\) must be known not only at 
the single point Xq but also in its vicinity (where Ta > 0, but small) in order to estimate 
the possibility of an experimental observation. 

Examples of T\(X) trajectories with T\(Xo) = as well as with Ta(Xq) w are studied 
on the basis of the FPO method for concrete systems, see [ljj] for atoms and 31 ] for quantum 



dots. In both cases, the trajectories depend strongly on energy near X and Xq, respectively. 
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This energy dependence is caused by widths bifurcation, i.e. by the avoided level crossing 
phenomenon. The interaction of the resonance states via the continuum (described by the 
complex non-diagonal matrix elements of the second term of H e g) plays an important role. 
The interplay between their real and imaginary parts makes possible the appearance of BICs 
at finite (physical) values of the coupling strength Every BIC appears together with 

at least one other state whose width is enhanced (due to the widths bifurcation) around 
X = X and whose energy is, generally, close to E\ . In the considered cases, the condition 
for the exact appearance of a (T symmetric) BIC is space reflection symmetry (V symmetry) 
of the system. Violation of V symmetry leads to T\(Xq) small but different from zero, i.e. 
to violation of time reflection symmetry (violation of T symmetry). 

The BICs, being VT symmetric states, are the result of widths bifurcation, i.e. finally of 
the existence of branch points in the complex plane. The branch points appear as families 
of trajectories when considered as a function of a parameter, as shown in numerical studies 
for double quantum dots 10]. As a consequence, also the BICs appear as families of tra- 
jectories by controlling the system by means of parameters. Besides the trajectories of the 
VT symmetric BICs there are always trajectories of states whose widths are large (widths 
bifurcation). These states are V symmetric, however T symmetry is broken as discussed in 
Sects. 3 and 4. 

It is interesting to remark that, using the condition of VT symmetry, Bender et al. 
obtained classes of non-Hermitian Hamilton operators whose spectra are real and positive 
. Recently, these solutions and their consequences for physical processes are discussed 
intensively in the literature, see e.g. the 2007 Workshop on Pseudo Hermitian Hamiltonians 



in Quantum Physics [351 ]. 



VI. SPECTROSCOPIC REORDERING PHENOMENA IN THE OVERLAPPING 
REGIME 

An unsolved problem in standard quantum mechanics is the description of the crossover 
from the regime with weak coupling between discrete and continuous states to that with 
strong coupling between them. For example, an interpolation procedure between the limiting 
cases with isolated resonances at low level density and narrow resonances at high level 



density (superposed by a smooth background term) is introduced in [37]] for the transmission 
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through a quantum dot. In contrast to such an interpolation procedure, the crossover can 
be described in the FPO formalism by means of the phase rigidity r\ that is reduced for 
many states A in this regime. 

Let us consider the one-channel case, C — 1, and ^ E —>■ ty E in the interior of the system. 
The right and left wave functions follow from (ITS]) and (1141) with df = cf * when excitation 
and decay of the state A occur via the same mechanism. Therefore the ^ E can be normalized, 

(^' left l^' right > = E c f 4 (<Wv> = El c fl 2 = 1 - (27) 

AA' A 

The normalization has to be done separately at every energy E due to the explicit energy 
dependence of the cf. Moreover, 

(^1^ = £cfcf, (0 A |0 V ) = E(cf) 2 A x (28) 

AA' A 

due to B% = -B%,, see ©. A x is a real number j4j. From ([27]) and ((2HJ) follows 

<*!l*!) =E(cf) ^ A = 5: M)!, (29) 

where the definition (ITU]) for r A is used. Then the phase rigidity p of the wavefunctions ^ E 
may be defined by 

P = ^ E = e™ E V (iMcf )] 2 - [Im(cf )] 2 ) (30) 

A ^ A ^ 7 

in analogy to (TTO]) . The value p corresponds to a rotation of ^ by # corresponding to the 
ratio between its real and imaginary parts. In spite of the complicated structure of p, it 
holds 1 > p > 0, see [8j]. The value p is uniquely determined by the spectroscopic properties 
of the system that are expressed by the coupling coefficients to the environment and the 
level density, or by the positions and widths of the resonance states and the phase rigidities 
30]) shows the relation between p and the r\. 



r x . Ec^ 

in y,y, the amplitude of the transmission through a quantum dot is considered in the 
framework of the S matrix theory, 

t = -2vuE ^ • ( 31 ) 

Here, the eigenvalues z x and eigenfunctions <p\ of H c g are involved with their full energy 
dependence, see Sect. 2. 
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For p = 1 and well isolated resonance states, the transmission amplitude (13T!) repeats 
the resonance structure of (fTTj) of the wave function ^/q. The transmission peaks appear 
at the positions E\ = Re(zx)\E=E x ~ of the resonance states. An analogous result 
holds when there is a nonvanishing background term additional to the resonance term f[3"Tj) 
of the transmission amplitude. The time scale corresponding to this so-called direct part 
of the transmission is, generally, well separated from that corresponding to the resonance 
part described by fl3T|) [3^]. Mostly, the resonances are narrow and well separated from 



one another. They appear as Fano resonances |40j on the smooth background. Due to the 
different time scales of the resonance and direct processes, it is \p\ ~ 1 also in this case. 

The situation is another one when the resonances overlap. In the overlapping regime, the 
resonance states avoid crossings with neighbored resonance states. In this case 

r A <4vr (iE\V\(l>x)mV\^) (32) 

holds even when there is only one channel in each of the two identical leads. The relation 
fl32|) differs from 

r A = 4vr ($\V\<f> x )(<f>* x \V\$) (33) 



' the eigenf unctions <p\ that can not be 
4j. Therefore, the contribution of the 



for isolated resonances due to the biorthogonality o: 
neglected in the regime of overlapping resonances 
state A to t^E^E x ) is larger than 1. The unitarity condition will be fulfilled, nevertheless, 
due to interferences and the possibility to rotate the <f)\, i.e. due to the non-rigidity of the 
phases of the wave functions <p\. As a consequence, the transmission in the overlapping 
regime does not show a resonance structure. Instead, it might be nearly plateau-like, for 



details see 



Bis, 



251 ] . Let us rewrite therefore the transmission amplitude fl3TT) by means of 



the wave function (TTBlh 

t = -2m (34) 

with being complex, in general. The advantage of this representation consists in the fact 
that it does not suggest the existence of resonance peaks in the transmission probability. 
Quite the contrary, the transmission is determined by the degree of alignment of the wave 
function ^ with the propagating modes in the leads, i.e. by the phase rigidity p. 
Nevertheless, the expressions fl34|) and fl3~T|) are fully equivalent. 
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The plateau-like structure of the transmission can not be obtained in standard quantum 
mechanics with fixed phases of the wave functions, r\ — 1 and p — 1. It is generated by 
interference processes with account of the alignment of some of the resonance states to the 
scattering states £q of the environment. At most, Re^^ = ±Im^^ (in the same manner as 
for the This case corresponds to p = 0. It will be reached when many resonance states 
are almost aligned with the £p, and J2\^ e [( c \E) 2 ]/r\ ~ according to (130]) . 



The numerical results 
method 



251 ] obtained by using the tight-binding lattice Green function 



361 ] for the transmission through microwave cavities of different shape show exactly 



the features discussed above. In the weak-coupling regime as well as in the strong-coupling 
regime, the transmission shows a resonance structure as expected from the standard quan- 
tum mechanics. The only difference between the two cases is the appearance of a smooth 
background term in the strong-coupling regime which does not exist in the weak-coupling 
case, and the reduction of the number of resonance peaks by two (corresponding to the align- 
ment of two resonance states each with one channel in each of the two identical attached 
leads) . 

In the crossover from the weak-coupling regime to the strong-coupling one, however, the 
calculated transmission is plateau-like instead of showing a resonance structure [251 ] . It 
is enhanced as compared to the transmission probability in the two borderline cases. In 
this regime, the resonance states overlap and spectroscopic reordering processes take place. 
Due to widths bifurcation, some of the resonance states become short-lived while other 
ones become trapped (long-lived). The enhancement of the transmission is caused by the 
short-lived states. Most interesting is the anticorrelation between transmission \t and phase 
rigidity \p\ which can be seen very clearly in all the numerical results obtained in 9|, |25J. The 
behavior of the transmission in the crossover regime with overlapping resonance states does 
not correspond to the expectations of the standard quantum mechanics with rigid phases of 
the eigenf unctions of a Hermitian Hamilton operator, and decay widths obtained from poles 
of the S matrix, see e.g. Moreover, the transmission in the crossover regime is not only 
enhanced but it also outspeeds the transmission calculated in standard quantum mechanics. 
The reason is the formation of aligned (short-lived) resonance states in the vicinity of branch 
points in the complex plane. 

The quantum brachistochrone problem of a physical system can be studied by considering 
the time needed for the transmission through the system from one of the attached leads to 
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another one |8(]. The transmission time at a certain energy E is determined by the delay 
time, i.e. by the lifetime of the resonance states lying at this energy. The lifetime of 
a resonance state is bounded from below: it can not be smaller than allowing traveling 
through the system in accordance with traveling through the attached leads, i.e. the system 
may become transparent at most [41]. This lower bound can be reached in a system with 
the non-Hermitian Hamilton operator H c g by aligning the wave functions of the system with 
those of the environment while such a possibility does not exist when H e g is considered to 
be Hermitian. 



VII. CONCLUSIONS 



In the present paper it is shown that the FPO technique is a powerful method for the 
description of open quantum systems. The core of the method is the definition of two 
subspaces and the non-Hermitian symmetric Hamilton operator H e g describing the localized 
states in one of the two subspaces. It guarantees a unified description of, on the one hand, 
discrete and resonance states and, on the other hand, resonance and decay processes. The 
spectroscopic information is received from the eigenvalues z\ and eigenfunctions <f)\ of H e R. 
Generally the eigenvalues are complex (with the exception of BICs and low-lying discrete 
states). In order to receive the spectroscopic information, the two subspaces must be defined 
in a meaningful manner. Observable values that are related to the S matrix and to the wave 
functions ^/^ (being solutions of the Schrodinger equation with the Hermitian Hamilton 
operator H) are independent of the definition of the two subspaces as long as P + Q = 1 is 
fulfilled. The S matrix is always unitary, and it is not necessary to consider its poles. 

The branch points in the complex plane determine the trajectories of the eigenvalues 
z\. level repulsion and widths bifurcation, respectively, appear in approaching them under 
different conditions. This phenomenon (avoided level crossing in the complex plane) causes 
the appearance of BICs in V symmetric systems, i.e. of VT symmetric solutions of the 
Schrodinger equation with non-Hermitian Hamilton operator. The BICs are localized at all 
times although no selection rule forbids their decay. Furthermore, the phases of the (f)\ are 
not rigid in the vicinity of the branch points. Due to this phenomenon, the eigenfunctions 
of H e Q may align to the scattering wave functions of the environment with the consequence 
that observable values are enhanced in the regime of overlapping resonances (where many 
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branch points exist). 

The branch points in the complex plane influence also the properties of discrete states: 
the phenomenon of avoided crossing of discrete levels known for a long time, can be traced 
back to the branch points. Generally, the branch points in the complex plane introduce 
nonlinear effects into quantum mechanics. This fact allows, among others, to understand 
the well-known relation between avoided level crossings and quantum chaos. 

Most interesting is that, due to the non-Hermiticity of H e g, some problems of standard 
quantum mechanics are solved without any additional assumptions. For example, time 
asymmetry appears in a natural manner since if e ff contains a time operator. The results 
obtained by using the FPO formalism with the non-Hermitian Hamilton operator H s r may 
be important for applications as well. By controlling the interplay between internal and 
external mixing of the resonance states by means of external parameters, the position of 
the localized BICs can be varied. Also the transmission through a system can be controlled 
parametrically: a system may become transparent in the crossover from the weak-coupling 
regime to the strong-coupling one. 
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